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Exact and quasi-exact solvabilities of the one-dimensional Schrodinger equation are discussed 
from a unified viewpoint based on the prepotential together with Bethe ansatz equations. This is a 
constructive approach which gives the potential as well as the eigenfunctions and eigenvalues simul- 
taneously. The novel feature of the present work is the realization that both exact and quasi-exact 
solvabilities can be solely classified by two integers, the degrees of two polynomials which determine 
the change of variable and the zero-th order prepotential. Most of the well-known exactly and 
quasi-exactly solvable models, and many new quasi-exactly solvable ones, can be generated by ap- 
propriately choosing the two polynomials. This approach can be easily extended to the constructions 
^^ ' of exactly and quasi-exactly solvable Dirac, Pauli, and Fokker-Planck equations. 
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^H ■ I. INTRODUCTION 

>0 Two decades ago, a new class of potentials which are intermediate to exactly solvable (ES) potentials and non 

solvable ones have been found for the Schrodinger equation. These are called quasi-exactly solvable (QES) models 
,__, ' for which it is possible to determine algebraically a part of the spectrum but not the whole spectrum [l|, 0, S, 0, [H, @, 
0, H B S3, El LL2, uM ■ The discovery of this class of spectral problems has greatly enlarged the number of physical 
Q-r systems which we can study analytically. In the last few years, QES theory has also been extended to the Pauli [14| 
and Dirac equations [HI, UM, LL3, [la LLiA L20, l22i [22, [231 . More recently, we have considered QES quasinormal modes, 
which are damping modes with complex eigen-energies |24J . Such modes are of interest in black hole physics. 

Usually a QES problem admits a certain underlying Lie algebraic symmetry which is responsible for the quasi-exact 

solutions. Such underlying symmetry is most easily studied in the Lie-algebraic approach [3. [a. M. la. M. flfJ. [ill Ha fl3| . 

However, solutions of QES states are more directly found in the analytic approach based on the Bethe ansatz equations 

rf\ . U» S [25| ■ In this analytic approach the form of the wave functions containing some parameters are assumed from the 

^. ■ very beginning, and these parameters are fitted to make the ansatz compatible with the potential under consideration 

Q\ \ ■ Further developments in QES theory include classification of one-dimensional QES operators possessing finitc- 

0"N dimensional invariant subspace with a basis of monomials [2b , and formulation extending to nonlinear operators 

^O ; [23. 

A different direction in the development of QES theory is the prepotential approach [28|] . Here the emphasis was 
t-H ■ shifted from the potential to the so-called prepotential (or superpotential) , a concept which plays a fundamental role 
in supersymmetric quantum mechani cs [29l. l3fj| (we mention here that the prepotential W considered here is the 
integral of the superpotetial W in [lj, [211 122, 1221, 1231 ) ■ Pre pot ential has been extensively employed to study classical 
and quantum integrability in Calogero-Moser systems [3l|, [32, [33, [34], (3f| [36| . The simplest prepotentials are those 
which give rise to QES potentials admitting just the ground states. Physically these are factors which take care of 
the asymptotic behaviors of the systems. They are called the gauge factors in [2, 3] (which we shall call the zero-th 
order prepotential). Classification of all possible gauge factors for sZ(2)-based QES systems is presented in Q- Unlike 
previous works, however, in the prepotential approach of 28] the prepotential assumes a more fundamental role. From 
a different consideration, it was found that the Schrodinger QES theory was most easily extended to systems with 
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multi-component wave functions by recasting the Lie-algebraic theory in terms of the prepotentials 1J, 



21, 

37] 



More 



recently, QES theory was extended to the Fokker-Planck equations also via the prepotential approach 

It thus seems that the prepotential approach to QES theory deserves a more in-depth study than it has received 
so far. The merit of this approach is that the form of the potential of the system concerned needs not be assumed 
from the beginning. All information about the system is contained in the prepotential and the solutions, or roots, 
of the Bethe ansatz equations. The prepotential and the roots determine the potential as well as the eigenfunctions 
and eigenvalues simultaneously. Also, in this approach exact and quasi-exact solvabilities can be treated on the same 
footing. Furthermore, as mentioned in the last paragraph, such approach facilitates extension of the QES theory from 
the Schrodinger equation to equations for multi-component wave functions. 

The ideas of the prepotential approach to exact and quasi-exact solvabilities have been presented in [l4|, [2l|, [22, (28| , 
and summarized in [33] . The emphasis of these works was placed on the feasibility and elegance of the prepotential 
approach. However, the forms of the the prepotential and the required change of coordinates were either directly 
adapted from the known ES and QES models, or given as known for the new QES systems. Now we would like to 
give a first attempt to address the questions as to how the choice of coordinate transformation and the prepotential 



are determined, at least for certain classes of coordinates and prepotentials. 

The main result of this work is the realization that exact and quasi-exact solvabilities can be solely classified by two 
integers, the degrees of the two polynomials which determine the change of variable and the zero-th order prepotential. 
Such classification scheme has not been explored before. There are upper limits for theses two integers beyond which 
no ES systems are possible. By selecting appropriate values of these two degrees, most of the well-known ES and QES 
models, and many new QES ones, are easily generated. This approach does not rely on the knowledge of possible 
existence of any underlying symmetry in the system concerned. It treats both exact and quasi-exact solvabilities on 
the same footing. 

This paper is organized as follows. In Sect. II we present the idea of the prepotential approach to the exact/quasi- 
exact solvability of the Schrodnger equation. The general conditions which the required change of variable and the 
choice of the prepotential are discussed. Two common types of transformations of coordinates leading to exact/quasi- 
exact solvable systems are discussed and examples presented in Sect. Ill and IV. In Sect. V and VI the prepotential 
is generalized to systems defined on the half-line and on finite interval, respectively. Sect. VII concludes the paper. 

II. PREPOTENTIAL APPROACH 

Suppose 4>o{x) (— oo < x < oo)) is the ground state, with zero energy, of a Hamiltonian H$: Hq4>q = 0. By the well 
known oscillation theorem <fio is nodeless, and thus can be written as 4>o = e~ w °^ x ' , where Wo (a;) is a regular function 
of x (this will be assumed in the rest of this paper). For the square- integrable 4>q, this is the simplest example of 
quasi-exact solvability. This implies that the potential Vb is completely determined by Wq: Vq = Wq — Wq, and 
consequently, the Hamiltonian is factorizable (we adopt the unit system in which ft and the mass m of the particle 
are such that ft = 2m = 1): 

This fact can be considered as the very base of the factorization method [38l l39( and of supersymmetric quantum 
mechanics [29|, yCJ • We shall call Wq (x) the zero-th order prepotential. 

Consider now a wave function <pN (N > 0) which is related to <fio of Hq by </>jv = 4>q4>n, where 



>N 



(z - zi)(z- z 2 )- ■ -{z- z N ), 00 = 1. (2) 



Here z = z(x) is some function of x. In taking the form of 4>n in Eq. ([2]) we have assumed that the only singularities 
of the system are z — ±oo. Other situations will be addressed later. The function (f>N is a polynomial in an (N + 1)- 
dimensional Hilbert space with the basis (1, z, z 2 , . . . , z ). One can rewrite 4>n as 

4> N = exp(-W N (x,{z k })) , (3) 

with the N-th order prepotential Wn being defined by 

A' 

W N {x, {z k }) = W (x) -J2 ln \ z ( x ) - z k\- (4) 

fc=i 

Operating on 4>n by the operator ~d 2 /dx 2 results in a Schrodinger equation Hn4>n — 0, where 

H N = -^ + V N , (5) 

V N = Wfi-W'Jj. (6) 

From Eq. |gj), V N has the form V N = V + AV N , where V = W^ - Wq, and 

Here the prime denotes derivative w.r.t. the variable x. 

AVjv is generally a meromorphic function of z with at most simple poles. Let us demand that the residues of the 
simple poles, z k , k = 1, . . . , N should all vanish. This will result in a set of algebraic equations which the parameters 



{zk} must satisfy. These equations are called the Bethe ansatz equations for {zk}- With {zk} satisfying the Bcthe 
ansatz equations, AVJv will have no simple poles at {z^} but it still generally depends on {z^}. Thus the form of Vn 
is determined by the choice of z' 2 and Wq{x). 

In what follows we would like to demonstrate that the choice of z 12 and Wqz' determine the nature of solvability 
of the quantal system. We shall restrict our consideration only to those cases where WqZ 1 — P m (z), z 12 — Q n (z) and 
z" = Ri(z) are polynomials in z of degree m, n and I, respectively. Of course, Q n (z) and Ri{z) are not independent. 
In fact, from 2z" = dz /2 /dz we have I = n — 1 and 



, . , _ ldQ n (z) 



z " = R n _ l { z ) = -^H. (8) 

Consequently, the variables x and z are related by (we assume z{x) is invertible for practical purposes) 

f* dz 

x(z) = ± / , (9) 

J VQn(z) 

and the prepotential Wq(x) is determined as 

W (x)= [ X dx(-^$=) . (10) 



y/Qn{z) 



z{x) 



Eqs. (HJ) and (fT0| define the transformation z(x) and the corresponding prepotential Wq(x). Thus, P m (z) and Q n (z) 
determine the quantum system. Of course, the choice of P m and Q n must be such that Wq derived from (|10p must 
ensure normalizability of <po = exp(— Wo). 

Now depending on the degrees of the polynomials P m and Q n , we have the following situations: 

(i) if max{m, n— 1} < 1, then in Vn(x) the parameter N and the roots z^'s will only appear as an additive constant 
and not in any term involving powers of z. Such system is then exactly solvable; 

(ii) if max{m, n — 1} = 2, then N will appear in the first power term in z, but z/c's only in an additive term. This 
system then belongs to the so-called type 1 QES system defined in Q, i.e., for each N > 0, Vn admits N + 1 
solvable states with the eigenvalues being given by the N + 1 sets of roots Zk's. This is the main type of QES 
systems considered in the literature; 

(iii) if min{m, n — 1} > 3, then not only TV but also z^s will appear in terms involving powers of z. This means 
that for each N > 0, there are N + 1 different potentials Vn, differing in several parameters in terms involving 
powers of z, have the same eigenvalue (when the additive constant, or the zero point, is appropriately adjusted). 
When Zfc's appear only in the first power term in z, such systems are called type 2 QES systems in Q. We see 
that QES models of higher types are possible. 

We will illustrate these general situations with specific examples in the following sections. For definiteness, in this 
paper we will only consider cases with z' 2 = Qi{z) = q2Z 2 + q\z + qo (i.e., z" — q2Z + qi/2), where q2, qi and qo 
are real constants (by taking n = 2 here we mean to include Q\ and Qo as special cases if some of the coefficients 
vanish). This choice of z' 2 covers most of the known ES shape-invariant potentials in [29, 30] and the sZ(2)-based QES 
systems in Q , and a new one discussed in [28|, [37| • Such coordinates are called "sinusoidal coordinates" , which include 
quadratic polynomials, trigonometric, hyperbolic, and exponential types. The connection of sinusoidal coordinates 
with ES theory has been extensively discussed in |40j |. 

With this choice of z 12 , we have 
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In deriving Eq. (fTTj) use has been made of the following identities 

N „ N 



£ 7 V( T - 2 £ " J - ( ^V ) ■ ( 12 ) 

nfi \z-zk){z-zi) -— - - - - > •■ - -' ' 
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-^ (z- Zk)(z- zi) 
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V Zk - zi 
( z k 


" fa z ~ Zk 


\Zk~ Zi 



(13) 
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E? V T = 2 £^-(-^) +*(*-!). (14) 

z -^ i (z - z k )(z - zi) <— - - ~> \ ■' - ~> I 
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Demanding the residues at z k 's vanish gives the set of Bethe ansatz equations 

P m (z k )-%z k -^-J2 9 ^ hl = 0, k=l,2,...,N, (15) 

1 4 ijtk Zk Zl 

Putting back the set of roots z k into Eq. (jlip. we obtain a potential Vn(x) without simple poles. The degree of P m {z) 
determines the nature of the solvability of the system, namely, for m — 1, 2, 3, . . ., the system is, respectively, ES, type 
1 QES, type 2 QES, and higher types QES, as discussed generally in the last section. 

To proceed further, we must specify P m (z). We shall discuss cases where z' 2 is linear and quadratic in z separately. 

III. EXAMPLES WITH z' 2 (x) = Qi(z) 

Let us first consider the case where z' 2 (x) = Q\(z), i.e., q 2 — 0. For defmiteness we take Qx(z) = 4Az + q n 
(qi = 4A), where A, go are real constants. This implies z"(x) = Ro(z) = 2A. Hence, the general solution of z(x) is a 
quadratic form of x: 

z(x) = Ax 2 + Bx + C, A, B, C: real constants. (16) 

It is easily checked that qo is related to A, B and C through z' 2 — 4Az + B 2 — 4AC. 

We now illustrate how some ES and QES models can be constructed in the prepotential approach by taking different 
values of m. 

A. Exactly solvable cases: m — 1 

Suppose rn = 1 so that P\{z) = A\z + Ao (A\, Ao real). By writing Pi — A\(z — Zk) + A\z k + Aq, one obtains 
from Eq. (fTT]) that 



N N ( 



AV N = -24^1-2^— ) Pl ( Zk )- A -y2 — 
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m Zk ~ Zl 



= -2A X N. (17) 

The last term in braces in Eq. (|LT[) vanishes when Zk's satisfy the Bethe ansatz equations (fl~5|) . Now N only appears 
as a parameter in an additive term, and not in terms involving powers of z in AVjy. The roots z^s do not appear at 
all. The additive term can be treated as the eigenvalue. The Schrodinger equation reads 



( r ^ + K 2 -<)e- WN =2A 1 Ne- 



w N 



(18) 



We see that the potential Wq 2 — Wq is ES: by varying N, one obtains all the eigenvalues 2A\N and the eigenfunctions 
4> N = exp(-Wiv). 

As an example, let us take Qi(z) — 1 and P\{z) — bz. A particular solution of z is z(x) — x. From Eq. (|10p one 
gets 



f x bx bx 2 

Wq(x) — / —dx = — — h const. 



(19) 



The integration constant is to be determined by normalization of the wave function. For 0o = exp(— Wo) to be 
square-integrable, one must assume b > 0. The BAE (fTS"]) are: 

bx k -Y] — - — =0, k = l,...,N, (20) 

*-^ Xh — Xl 

ijtk K ' 



and the potential is 

V = b 2 x 2 - 6, AV N = -2Nb. (21) 

This leads to the Schrodinger equation: 

jL+b 2 xAe- WN = b(2N + l)e- WN . (22) 

ax 1 * / 

This system is just the well-known simple harmonic oscillator. 

We note here that by rescaling vbxk — > Xk, Eq. ([20)) will have 6=1. The resulted equations are the equations that 
determine the zeros of the Hermite polynomials Hn(x) as found by Stieltjes [3l|, [4j], [42] • Hence we have reproduced 
the well known wave functions for the harmonic oscillator, namely, 4>n = exp(— Wn) ~ exp(— bx 2 /2)H^(-Jbx). 

B. Type 1 quasi-exactly solvable cases: m — 2 

Next we consider P2(z) — A2Z 2 + A\z + Aq. By a similar argument we obtain 

N 

AV N = -2A 2 Nz - 2A 2 ^ z k - 2A 1 N. (23) 

k=l 

The Schrodinger can be written as 

~ + W 12 ~ W " - 2A 2 Nz\ e- WN = 2 I A 2 J2*k + A^) e~ WN . (24) 

Unlike the previous case, now TV not only appears in an additive constant term but also in the term with z, and the 
set of roots Zfc's appear in the additive term. This system is the so-called type 1 QES models. Type 1 QES models 
classified as class VI in [3[ belong to this category. 

A well-known example is the sextic oscillator, the simplest QES model of this type [l|. In our prepotential approach, 
this system is defined by z(x) — x 2 and P 2 (z) = 2(az 2 + bz). Then Q\(z) — 4z, and 

= r x ax +bx ^ = l^ 1^2 congt 
J Vx 2 4 2 

Here a > to ensure square-integrability of the wave function. The BAE are: 

2azl + 2bz k - 1 - 4Y^ Zk =0, k=l,...,N, (26) 

i^k Zk ~ Zl 

and the potential is 

V N = a 2 x 6 + 2abx 4 + [b 2 - (AN + 3) a] x 2 - Aa ^ z k - (AN + 1)6. (27) 

fc 

It is seen that the QES sextic oscillator can be so easily constructed in the prepotential approach. 

C. Type 2 quasi-exactly solvable cases: m = 3 

We now consider cases with m — 3 with P%(z) — A3Z 3 + A 2 z 2 + A\z + Aq. Eq. (fTT|) leads to 

(N \ N N 

A 3 J2 z k + A 2 N\z- 2A 3 J24- 4A 2 J2 z k ~ 2A X N. (28) 

fe=i / fc=i fc=i 

Now N appears in z and z 2 terms, and also in an additive constant term. The roots z^'s now not only appear in the 
additive term but also in the term with z. This is a type 2 QES model. 



A simple example of this type is given by the defining relations z(x) = x and Pz{z) = az 3 + bz. The prepotential 
is Wq — ax A /4 + bx 2 /2, which is exactly the same as that for the sextic oscillator discussed in the last section. The 
two models differ only in the choice of z(x). The corresponding Vjv is 



N \ N 



V N = a 2 x 6 + 2abx A + [b 2 - (2N + 3) a] x 2 - 2a \Y^ x k \x - 2a ^ x\ - (2N + 1)6. (29) 

\fc=i / fc=i 

This is a new QES model. 

It is now easy to see that, for m > 4, not only will A appear in more terms involving powers of z, but also the set 
of roots Zk's. This will give rise to new general types of QES systems as mentioned in Sect. II. 

IV. EXAMPLES WITH z' 2 (x) = Q 2 (z) 

We now come to cases where z' 2 (x) = Qziz) with q2 7^ 0. Here z(x) is again some sinusoidal coordinates, which 
include the exponential, hyperbolic and trigonometric functions. Construction of models proceeds as before. By 
taking appropriate P m (z), one can reconstruct class I. II and X QES models in 3], and some of the ES models listed 
in [23|, namely, the Morse potential, the Scarf I and II potentials, and the Poschl- Teller potential. We will not bore 
the reader by going through all the cases here. Instead we shall briefly discuss the Morse potential, as we would like 
to show how easy its potential, eigenvalues and eigenfunctions are constructed, and to compare this construction with 
another construction based on a different choice of the prepotential to be discussed in the next section. 

Suppose we take Qi{z) = a 2 z 2 , and choose a solution z(x) — exp(ax) (henceforth a is taken as a positive real 
constant). Let P\(z) = a(Az — B) (A, B > real constants). The parametrization is chosen such that the form of 
the Morse potential given in [29| is recovered. From Eq. (fTO]) we have Wq = A — B exp(-ax). Hence 

V N =A 2 -B (2A + a) er ax + B 2 e - 2ax - \a 2 -(A- Na) 2 } , (30) 

with Zfc's satisfying the BAE 

a v— v z? 

Az k -B--z k -aY] *— =0. (31) 

This is the ES shape- invariant Morse potential listed in [29| . Taking the first three terms in Eq. (f3"0"|) as the traditional 
Morse potential, the eigenvalues are given by A 2 — (A — Na) 2 , in agreement with the result given in[29]. The wave 
functions are 

/ R \ N 

b N (x)~exp[-Ax--e- ax \ ]J(z - z k ), z = e ax . (32) 

' fe=i 



n 



Let us make an interesting observation here. From Eq. (TTTJ) it is obvious that only the coefficients qi in Q<z{z) 
and pi in P\(z) will enter the expression of the eigenvalues, namely, N(2p\ ~ q^N). Now for the Scarf II and the 
Poschl- Teller potential (both with Qz(z) = a 2 (l + z 2 )), qi = a 2 is the same as that in the Morse potential, whereas 
for the Scarf I potential (with Q<2,{z) — a 2 (l — z 2 )) there is a sign difference. So if we choose P\{z) with the same p\ 
for all these cases, e.g. P\{z) = aAz +po, then we would expect that the Morse, Scarf II and Poschl- Teller potentials 
would have the same set of eigenvalues (A 2 — (A— Na) 2 ), while the Scarf I potential has a different set of eigenvalues 
differing by a sign in some parameter ((A + Na) 2 — A 2 ). This is in fact the case [29J. The prepotential approach 
presented here gives a very simple and direct explanation of why this is so. 

V. PREPOTENTIALS FOR SYSTEMS WITH SINGULARITY z = 

Now we would like to discuss a possible generalization of the A-th order prepotential in Eq. Q for quantum systems 
defined on a half-line (e.g. x £ (0,oo)) with singularity at the origin. For such systems, the wave functions may 
acquire a prefactor x p , where p is usually some non- negative positive number, in order to account for the asymptotic 
behavior at the origin. This observation motivates a possible generalization of the zero-th order prepotential Wq(x) 



to Wq(x) = Wq(x) — p\n\z\ for the ground state 4>q(x) = exp(— Wo(x)). Here Wo is a regular function of a; as before. 
Eq. ([4]) becomes 



N 



W N {x,{z k })=W (x)-phi\z\-^2ln\z(x)-z k \, z k £ 0. (33) 

fe=i 

For the moment p is a free parameter. 

With the prepotential (|33|) . the potential Vat = Wjy — W'^ has the form V~o + AVn where Vo and A Vat are given by 

V = W' 2 -W' a ' 

= <-<- 2 (^'-^)|+p(p-l)(^ 2 (34) 

AViv = -2(w>z'-^)y2^— + z> 2 V 7 1 - + V , 2P A . (35) 

1 '' ' ^ - ■ ' X? ( z - z k){z - zi) f^ziz-zk)} 



2 J ■f—'z - z k 

k—l 



\> 



Again, the system is completely defined by the choice of z' 2 and Wo- However, the presence of terms with the 
parameter p changes qualitatively the sufficient conditions discussed in Sect. II and opens up new possibilities. 

As in previous sections, for definiteness, we shall confine our discussions here to z' 2 = Q%(z) = q2Z 2 + q\z + go and 
Wqz' = P m (z). With these choices, the potential Vat is given by 

Vat = W^ 2 ~W^-p(2P m {z)- qi (p-]^\- z+ p{p-l) q() ^ + q2 (N+p) 2 




H)-g£*S- ^ 



The corresponding Bethe ansatz equations are 



PM + ^ - (p+ I) q ,z k -( p+ l) qi -J2^=0. (37) 



z k \ 2j v 4; i*k Zk - zi 

It is clear that Eqs. ||3BJ) and $57} reduce to (dJ) and JT5J when p = 0. 

If go = (i.e., z = is a zero of z' 2 ) and m < 1 one may obtain an ES model. For example, if we take Qi(z) = az 
and P\{z) = A0, then the three-dimensional oscillator listed in [29J is recovered. But when <?o ^ 0, the presence of 
the pqo/ z{z ~ z k ) term in Eq. (|36p will give rise to a term (J^ fe Zk)/z in Vjv, and hence the potential (J3SJ) represents a 
QES system even if m < 1, in contrast to the cases discussed in previous sections. For instance, if we take z' 2 = 1 and 
Pm{z) — P\(z) = a+bz, this will produce an ES model if p = 0, which is the shifted oscillator Vjv = (bx+a) 2 — (2N+l)b 
(note that the domain of x changes from the half-line to the full line) . But for general p, the system is a type 2 QES 
system classified as the class VIII system in [3j. 

However, it is possible that if p assumes certain value, the nature of the system could be qualitatively changed, 
such as the domain of the variable may change from the half-line to the full line (as in the case of the shifted oscillator 
mention in the previous paragraph), or a QES system becomes an ES one. We shall illustrate these situations with 
two examples below. 

A. Sextic oscillator again 

Following Sect. III(B), we take Q<2.{z) — 4z and P m (z) — 2(az 2 + bz). These lead to z — x 2 and W (x) — 
aa; 4 /4 + bx 2 /2 + c, with real constants a > 0, b and c. The BAE and Vn are 

2azl + 2bz k - Up + 1) - 4 V Zk = 0, k = l,...,N, (38) 

L — ' Zk — zi 
lj£k 

and 

V N = a 2 x 6 + 2abx 4 + [b 2 - {AN + Ap + 3) a] x 2 

+Ap(p-^\^~AaY / Zk-^N + Ap+l)b. (39) 



In general, this model is a QES system on the half-line. But if p = or p = 1/2 the 1/x 2 term will be absent, and 
the system is just the sextic oscillator on the full line. The domain is extended. The case p = with symmetric wave 
functions was discussed before, which is the class VII QES model in J3|], and the case p = 1/2 with anti-symmetric 
wave functions was discussed in [l|, 0, [f| . 

B. Morse potential again 

Now we consider the situation which will include the ES Morse potential as a special case. 

We take Q2{z) — a 2 z 2 and P m (z) = Pi(z) = a 2 z 2 /2 — aAz. As a special case, we choose the solution z — exp(— ax) 
and Wq(x) = A — Bz with B = a/2. As in Sect. IV, the parametrization is chosen such that when p assumes special 
value, the system becomes the Morse potential given in [29|]. The potential is 

V N = A 2 - B (2A + a) e - ax + B 2 e~ 2ax - 2pa (-z - A\ + a 2 (N + pf 

If we choose z^ 's to satisfy the BAE (|57|) 

^ Zk -A-a(p + ^\-aJ2^— = Q> k = l,2,...,N, (41) 

^ ' l^k Zk Zl 

then we arrive at the potential 

V N = A 2 -B(2A + a)e- ax + B 2 e- 2ax -a 2 (N + p)z 

-a 2 ^2z k + 2aA{N + p) + a 2 {N+pf . (42) 

k 

The term — a 2 J^k Zk can ^ e simplified using Eq. (I41|) as 

] ^2z k = -a [2A + (2p + 1) a] N - a 2 N (N - 1) . (43) 



a 2 



In general this potential defines a type 1 QES system, as N appears in the first power term of z. This system is not 
listed in ||. But if p = 0, then the domain of x changes from the half- line to the full line, and the system becomes 
that classified as class I in |3|. 

It is also obvious that if p = —N, Vn becomes Eq. (13T)|) . and the system is the ES Morse potential. Let us recast 
Eq. (14"1"T) into the following form (taking p = —N): 

E^ + — = 1 k=l,2,...,N, (44) 



Zk — zi Zk 



where 



7=2 2N+1. (45) 



Eq. pi)l is just the set of equations that determines the zeros of the associated Laguerre polynomials Lj^ (z), i.e., 
L~jf (zk) = [3lll42|. Hence the eigenfunctions 4>n — exp(— Wn) are 

/ B \ N 

4> N (x) = exp(-Ax--e- ax )z- N Y[(z-z k ) 

^ a ' fe=i 

= z s - N e-Hf- N \z), a=-, z = e-°*, (46) 



JN 



as given in [2£ 



Compared with the discussion in Sect. IV, it is interesting to see that the Morse potential can be constructed with 
two different prepotentials. Now one may wonder if the results are consistent, as the wave functions and the BAE 
look rather different in these two constructions. Below we would like to show that they are indeed the same. 

Let us rewrite Eq. (j3"2"|) as 



4>n(x) ~ exp ( — Ax e 



N 



- ■ \ ~JV 



n^ -1 -** 1 ). z = 



(47) 



fc=i 



We note that in Sect. IV the variable z = exp(aa;) is reciprocal to the variable z = cxp(— ax) in this subsection. 
Hence when one makes the change z — > \/z and z k — > l/z k in Eq. (|47p . one arrives at Eq. (|46[) . Now one needs only 
to show that the same transformation in z and Zk maps the BAE (f5T|) to Eq. (|4T|) . 
Making the change Zk — > 1/ z k in Eq. (|3~lj) leads to 



A---Bz k -aJ2 — 

2 U zi Zk 



= 0. 



(48) 



Writing Z\ = (zi — Zk) + z k in the numerator of the last term in Eq. (|48[) and recalling that 5 = a/2, we arrive at 

1\ TT^ z k 



-z k -A-al-N +2 



E 



Zk -zi 



0, fc = l,2, . 



,JV. 



(49) 



This is simply Eq. (|41|1 with p = — JV. Thus we have shown that the wave functions and the BAE obtained in the two 
constructions are the same. 



VI. SYSTEMS DEFINED ON FINITE INTERVALS 



Finally, we consider systems defined on a finite interval. Suppose the potential of a system is singular at z = ai 
and 02, where a\ and a 2 are two real parameters (assuming a 2 > a\). Generalizing the discussions in the last section, 
it is plausible to assume for such system a prepotential of the form 



JV 

W N (x,{z k }) = W (x) - pi\n\z - ax\ - p 2 ln\z - a 2 \ -^Iti\z(x) - z k \, z k ^a\, a 2 , 

fe=i 

where p\ and p 2 are two real positive parameters. The wave function has the form 

JV 

<f> N ~ exp(-W (x))(z - a x ) Pl (z - a 2 ) P2 ]J(z z k ) 

Now Eqs. (f34|) and (|35|) are generalized to 

V = W^ 2 - Wtf - 2 (V^' - ^ / '"' 



fc=i 



P2 



+z- 



2 J \z — a\ z — a 2 
pi(pi-l) p 2 (p 2 -l) , 2pip 2 



(z-ai) 2 (z-a 2 ) 2 (z-ai)(z-a 2 ) 



AV/ 



iY 



'k-i)^ 



2 J f-j z - z/, 

fc— 1 



+z 



'2 



y * 



pi 



i>2 



z - ai z - a 2 j f— | z - z fe 



" 1 

Et 2 - 



(50) 



(51) 



(52) 



(53) 



Once again, the system is completely defined by the choice of z n and Wq. The analysis of the solvability of the system 
is more complicated than before, but the principle is the same. For polynomial z' 2 and Wqz' the systems are new 
QES models in general. 
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As before let us take z 12 to be at most quadratic in z. A situation of interest is that in which a\ and a 2 are the two 
real zeros of Q2(z), i.e., if z 12 = A(z — a-\){z — z 2 ) where A is real. In this case the sufficient conditions for ES and 
QES models are the same as before, namely, one could get ES system if m < 1, and QES otherwise. But just as the 
Morse potential discussed in the last section, the QES system could become an ES one if p\ and P2 take on special 
values. The ES Rosen-Morse I and II potentials and the Eckart potential in [29| are such cases. 

We shall illustrate the construction of a QES model first discussed in [28| (see also [331 )■ This model was not listed 
in (3[ but is found to be also related to sl(2) algebra (37|. The function z(x) is taken to be z(x) — sin 2 x. This is a 
solution of z 12 — 4z(l — z), and hence a\ = and a 2 — 1, and < x < ir/2. In order to obtain a type 1 QES model 
we should choose Wq such that Wqz' is of the second degree in z. Let us take P2(z) = 4az(z — 1), where a is real. 
This gives a solution Wo(x) = acos2x/2. With the chosen Q2{z) and P 2 {z), we obtain from Eq. (|B"2")) and ([53")) the 
following QES potential 

V N (x) = a 2 sin 2 2x + 2a cos 2x + 2 (a sin 2 2x + cos 2x) ( — ^ ^— ) 

\sin x cos^xj 

+4pi(>i - 1) cot 2 x + 4p 2 (p 2 - 1) tan 2 x (54) 

N 

-8aN sin 2 x - 8a^2 z k - AN [N - 2(a - pi -p 2 )] - 8p x p 2 , 
fc=i 

where the z^s satisfy the BAE (k = 1,2, ... ,N): 

Aaz 2 - 2 (2(o - Pi - p 2 ) - 1) z k - 1 - 4z fe (l - z k ) V ^^ = 0. (55) 

fe Zfe - z; 

The forms of V/v and BAE in [37| for this model is regained by identifying pi — c/2 — 6/4 and p 2 — b/4, where b and 
c are the parameters used in that paper. 

VII. SUMMARY 

We have discussed exact and quasi-exact solvability of the schrodinger equation based on the approach of prepo- 
tential. Three types of iV-th order prepotentials are described which cover most of the known ES and QES models, 
and which are capable of generating new ones. It is shown that the choice of coordinate transformation z' 2 {x) and 
the zero-th order prepotential Wo(x) completely determine the form of Vjy. General conditions on the choice of z' 2 
and Wq(x) for exact and quasi-exact solvabilities were described. These conditions were illustrated by various exam- 
ples. The prepotential approach is quite elegant in itself. Moreover, it allows easy extensions of ES and QES theory 
to systems with multi-component wave functions, such as the Pauli and Dirac equations, and to the Fokker-Planck 
equations, as prescribed in [14, I2U |22|. [37| . 

In this work we have confined our discussions to ES and QES models involving a change of coordinates to the 
so-called sinusoidal coordinates. These are coordinates z(x) which satisfy z' 2 = Qiiz), or z" = Rq(z) or R\(z). These 
coordinates cover most of the known ES and QES models. The examples we presented here by no means exhaust all 
possible sinusoidal coordinates. Other QES cases admitted by certain sinusoidal transformations need be studied. 

But as discussed in Sect. II, to construct a type 1 QES system it is also possible to choose z' 2 to be a cubic 
polynomial, i.e. z' 2 = Qs(z) . This lies beyond the transformations by sinusoidal coordinates. Such transformations 
are less studied in the literature. However, they could give rise to new QES systems, and deserve further study. In 
fact, two sZ(2)-based QES cases in [3j, namely, class IV and V, require n = 3 and m — 1 and 3, respectively. 

Only three types of JV-th order prepotentials were discussed here, namely, prepotentials for systems defined on the 
whole line, on the half-line, and on a finite interval. Generalizing to several singularities is straightforward. 

We have only considered cases for which Wqz' and z' 2 are polynomials in z. It is interesting to consider other 
possibilities for these two defining functions. 

Finally, it is also interesting to extend the present ap pro ach to systems with non-Hermitian Hamiltonians admitting 
real spectra. A preliminary attempt was presented in [43j, where some non-Hermitian QES Hamiltonians, including 
that given in [44j, were generated by making some coefficients in P m (z) complex. However . a full treatment of this 
kind of systems usually requires one to extend the basic variable x to the complex plane J4f|, [46[. Hence a better 
understanding of the prepotential approach to non-Hermitian Hamiltonians is needed. 
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